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We study a two-dimensional t — J model close to the Ising limit in which charge inhomogeneity 
is stabilized by an on-site potential Cs, by using diagonalization in a restricted Hilbert space and 
finite temperature Quantum Monte Carlo. Both site and bond centered stripes are considered and 
their similitudes and differences are analyzed. The amplitude of charge inhomogeneity is studied as 
Es — > 0. Moreover, we show that the anti-phase domain ordering occurs at a much lower temperature 
than the formation of charge inhomogeneities and charge localization. Hole-hole correlations indicate 
a metallic behavior of the stripes with no signs of hole attraction. Kinetic energies and current 
susceptibilities are computed and indications of charge localization are discussed. The study of the 
doping dependence in the range 0.083 < a; < 0.167 suggests that these features are characteristic of 
the whole underdoped region. 

PACS: 71.10.Fd, 74.80.-g, 74.72.-h 



I. INTRODUCTION 

After many years of intensive research the physics 
of underdoped high-Tc cuprates is still far from be- 
ing wholly understood. This underdoped region, where 
strongly correlated electron effects are most important, 
are considered by many to contain the key to explain 
central features of these materials, including the super- 
conducting phase. One of the main features of the 
underdoped cuprates is the presence of the pseudogap 
which has been-intensely studied with several experimen- 
tal techniques. EfH In spite of this intense effort, there is 
still not a clear picture of the origin of the pseudogap 
and its relation with the superconducting gap.u Another 
interesting phenomena observed by experiments in some 
underdoped cuprates is the presence of ordered charge 
inhomogeneities ( "stripes" ) followed (as the temperature 
is decreased) by the a]ppe|arance of incommensurate (IC) 
magnetic correlations These charge inhomogeneities 
provide a natural explanation for the pseudogap. There 
is also strong experimental evidence of the presence of 
charge ordering/stripes in otheE. iransition metal oxides 
like manganites and nickelates.Q'B Again, in spite of an 
enormous amount of research, there is not yet a clear 
picture of the mechanism leading to stripe formation and 
the relation between stripes and superconductivity in the 
cuprates. Different and in many cases competi]ig|-scenar- 
ios have been proposed to explain these issues. aii3 

From the theoretical point of view, in addition to 
the conceptual complexity of strongly correlated electron 
systems, there are important methodological difhculties. 
In stripes theories, phenomenological models of coupled. 
Luttinger (or Luther-Emery) liquids have been studiecO 
providing useful insight in the problem although quanti- 
tative estimations are usually lacking. In any case, these 
phenomenological models have to be supported by the 
study of models at a more microscopic level. These mi- 



croscopic models (variants of Hubbard, t — J models) are 
most frequently studied by numerical methods which also 
have to deal with severe limitations (some of them will 
be discussed in Section ||). It is then absolutely neces- 
sary to complement different approaches and to contrast 
the results obtained by analytical and numerical tech- 
niques in order to get a consistent picture of the physics 
of underdoped cuprates. 

In this article a model of stripes is studied by numerical 
techniques. The model is defined as a two-dimensional 
(2D) t — J model where stripes are stabilized by an effec- 
tive on-site potential representing the effect of extrinsic 
features to Cu-0 planes like lattice anisotropy,Ej effect 
of chains in YBa2Cu3 07_i (Rcf. |^) or other structural 
details, Coulomb interactian due to off-plane charges, 
electron-phonon coupling,t3 etc. This model has been 
stud ie | ti .numerically on small clusters (4 x 4) at both 
zerotSliB and finite temperatures .liH. Results on larger 
clusters with periodic boundary conditions were obtained 
in Ref. although at the price of working close to the 
Ising limit of the Heisenberg term of the Hamiltonian. 
In the present article we present an exhaustive analysis 
of this model by using quantum Monte Carlo (QMC) 
techniques, in particular, we examine the behavior at 
various doping and the effects of the anisotropy of the 
exchange term of the Hamiltonian. These analysis at 
finite temperature are complemented with zero tempera- 
ture studies with a diagonalization technique in a reduced 
Hilbert space kimwn as "systematically expanded Hilbert 
space" (SEHS).N The resuffs obtained with both QMC 
and SEHS allow a comparison not only with experiment 
but also with previous calculations on similar rlusters but 
with (partially) open boundary conditions.ElL^ 

The paper is organized as follows. In Section ||, the 
model considered is defined and some elementary details 
of the techniques employed are discussed. In Section |l| 
we present results obtained with QMC and SEHS for the 
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stripe filling as a function of the on-site potential both for 
site- and bond-centered stripes (defined in Section ||) . In 
Section [V we study charge and spin ordering as a func- 
tion of temperature and in Section ^ we show the behav- 
ior of the kinetic energy and current susceptibility with 
temperature as an indication of charge transport prop- 
erties in the underdoped region. Finally, in the Conclu- 
sions we summarize the results as an unified picture of 
the evolution of the various physical quantities examined 
as a function of temperature. 



II. MODEL AND METHODS. 

The Hamiltonian of the t~ J model with an anisotropic 
exchange term is: 



Htj= - i 51 (^I<TCi<T + h.c.) 



+ Jll(l^^tSj + SrS+) + S:S^ - \n,n,) (1) 

where the notation is standard. i,j are sites on a square 
lattice. N site clusters with fully periodic boundary con- 
ditions (PBC) were considered. We adopted i = 1 as the 
unit of energy and temperature and J — 0.35. In some 
cases we considered J = 0.7, specially because it reduces 
the sign problem. 7 = corresponds to the Ising limit or 
t — Jz model.Ell The stripes are stabilized by an effective 
on-site potential: 



(2) 



On "site centered" (SC) stripes Csi — — < on equally 
spaced vertical columns, as in the original picture of 
Ref. ^, and e^j = on the remainder sites. At each 
doping fraction, the number of stripes imposed is such 
that the linear hole density on the stripes at infinite po- 
tential is 1/2. We shall come to this point later on this 
section. 

The QMC method employed is an implementation of 
the usual wocidpline algorithm with the checkethoard 
decompositioriE3'E3. The "minus sign" problemEj be- 
comes more severe as lattice size and doping are increased 
and temperature is decreased. It is absent if holes are 
confined to move in one dimension and the XY tern in 
the Heisenberg interaction is set to zero. The behav- 
ior of the average sign as a function of various param- 
eters of the model is shown in Fig. 1 of Ref. |l8|. We 
checked our code against for the 4x4 cluster with one 
site-centered stripe, 2 holes, J/t = 0.35, Cs — 1. We ex- 
trapolated our results to zero temperature and compared 
them against the ones obtained with Lanczos exact di- 
agonalization. The relative error for the energy was of 
the order of 10~^, and for correlations 10~^. For larger 
clusters, smaller e^, or larger 7 the correlation's error es- 
timated from different independent runs may be as large 



as 10 % at the lowest temperature attainable. In gen- 
eral, we kept A = (3/M = 0.083, where (3 = l/T and M 
is the Trotter number. For the lowest temperatures of 
our simulations {T/t « 0.1) the measurement runs were 
as long as 1.6 10^ MC sweeps. With these very long sim- 
ulations we were able to deal with {sgn) as small as 0.05. 
On the largest clusters we considered, 8x8 and 12 x 12, 
we studied ~ 2, which corresponds to well-defined 
charge inhomogeneities (see next sectimiL ^-^d has been 
also taken in previous related studiesO'EEl 

SEHS is a zero temperature diagonalization in a 
Hilbert space that is expanded by successive hole hop- 
pings, starting from an initial state (e.g. holes on the 
stripes) . After each expansion the Hamiltonian is diago- 
nalized with the Lanczos algorithm and a fraction of the 
most weighted configurations of the ground state is kept 
for the next iteration.EElThe goal is that after several back 
and forth steps, an optimal set of states is selected for a 
given dimension of the Hilbert space. The Hilbert space 
is grown to a maximum dimension (determined of course 
by the computational resources available) and finally the 
quantities of interest are extrapolated to the full dimen- 
sion of the Hilbert space of the system. In many cases 
reliable qualitative behaviors are obtained without this 
final extrapolation. The whole procedure is variational 
but it shares many of the virtues of exact diagonaliza- 
tion like the possibility of computing any physical quan- 
tity both static and dynamical. Another advantage that 
SEHS shares with Lanczos diagonalization is the possi- 
bility of working for a given set of quantum numbers, 
in particular since periodic boundary conditions can be 
used, it is possible to work at a given momentum. In 
addition there are no constraints on the model studied. 
In the following section we will apply this technique to 
the t — Jz model but in principle it could be applied to 
the isotropic exchange although in this case requirements 
of time and memory are much larger. Another possibil- 
ity offered by this technique is that long-range Coulomb 
interactions can be handled, as a difference with QMC. 
We will explore in the future this possibility. 

The use of the on-site potential to stabilize stripes 
could be thought as leading to "static" stripes, not in 
the sense of "pinned" (since we are working with peri- 
odic boundary conditions) but in the sense of suppress- 
ing stripe ffuctuations. The opposite meaning of "static" 
could be applied to the stripes resulting from dena-ti^ ma- 
trix renormalization group (DMRG) calculationst2l with 
open boundary conditions along the longer direction of 
the cluster. In this case, it is well known that "bond- 
centered" (BC) stripes are obtained rather than SC ones. 
This fact may refiect the also well known tendency of 
holes to forp, pairs in the 2D t — J model at zero 
temperatureEil However, SC stripes are much sharper 
features than BC stripes and they seem more corisistent 
with elastic neutron scattering, NQR and ARPESE3 data 
at low doping although not definite conclusions can be 
drawn yet. At higher doping there have been specula- 
tions that more BC stripes may be formed at the ex- 
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jC stripes thus explaining some experimental 



III. STRIPE FILLING. 

Our first study concerns the evolution of charge inho- 
mogeneity as a function of the applied on-site potential 
in the 2D t ~ Jz model. Zero temperature results were 
obtained with SEHS for a 8 x 4 cluster, 4 holes (hole 
doping fraction x ~ 1/8), PBC. Stripes were imposed 
along the shorter direction {y axis). Translations along 
both directions, reflection along an horizontal axis and 
spin reversal symmetries were imposed and the dimen- 
sion of the Hilbert space was grown up to w 8.6 10®. 
The ground state in the case of SC stripes corresponds 
to {qx,qy) — (0,7r) and to (0,0) in the BC case. 
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FIG. 1. Hole density on stripe sites on the 8x4 cluster, 
J/t = 0.35, as a function of the on-site potential e^. Circles 
(triangles) correspond to SC (BC) stripes. Solid and dashed 
lines are quadratic extrapolations to H=j 0. The star in- 
dicates the result obtained with DMRCEZl for (almost) BC 
stripes on the same system but mixed boundary conditions. 
The inset shows the stripe densities as a function of the di- 
mension of the Hilbert space and their extrapolation to the 
full dimension. From bottom to top = 0.2, . . . , 0.6. 

Results for the stripe filling are depicted in Fig. ^ for 
J/t — 0.35 as a function of the on-site potential. Details 
of the extrapolation to the full Hilbert space dimension 
for each eg with a power law are shown in the inset of 
this figure. These densities are then extrapolated with a 
quadratic polynomial to Cs = 0. The final extrapolated 
values of the stripe densities are very similar (0.153 for 
SC and 0.159 for BC stripes) and they are quite close 
to those obtained with DMRG in Refs. 0|2O[ Of course 



the linear density in the BC case would correspond to 
the double of that quantity, but this concept is some- 
what misleading because it hides its relative smoothness 



with respect to SC stripes. It has been also previously 
shown that-the effect of exchange anisotropy is not very 
important. E2l We also include for comparison the DMRG 
result for the same model, cluster and parameters but 
using mixed boundary conditions (open (periodic) along 
the direction perpendicular (parallel) to the stripes). lZI 
In this case, the occupancy of the two legs on a stripe is 
slightly different and their average value was taken. As 
expected, the open boundary conditions along the direc- 
tion transversal to the stripes leads to a somewhat larger 
charge inhomogeneity than our results with fully PBC. 
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FIG. 2. Hole density on stripe sites on the 8x4 and 8x8 
clusters, J/t = 0.35, as a function of the on-site potential 
es favoring a SC stripe. The dashed line corresponds to the 
uniform doping x = 1/8. Points above this line correspond to 
the stripe; below it, they correspond to the first column next 
to the stripe. Solid lines are quadratic fitting to the data. 

In the case of SC stripes we can compare the previous 
results with the ones obtained with QMC. The tempet^ 
ature dependence of the stripe density is very smoothed 
and it has been extrapolated to zero temperature with a 
power law. The stripe hole density on the 8x4 cluster 
obtained with QMC is shown in Fig. |[ The resulting 
error of the finite temperature data and of the extrapo- 
lation procedure can be appreciated from the dispersion 
of the data. Nevertheless, the overall agreement between 
QMC results and the ones obtained with SEHS is cer- 
tainly very good. The results for a larger cluster, 8x8, 
are quite consistent with the previous ones. The stripe 
density seems somewhat smaller than for 8 x 4, which 
could be explained by the geometric anisotropy of this 
cluster making it more "one-dimensional" thus favoring 
the formation of a charge density wave along the long 
axis. However, error bars can likely mask this effect, if it 
exists. 

From the results of Fig. || one can also infer that to 
have a stripe filling close to 1/2, let us say, larger than 
0.4, the effective on-site potential should be w 1 in the 
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case of SC stripes, and somewhat smaller (« 0.35) in 
the case of BC stripes. Then, a sizable on-site poten- 
tial, specially in the case of SC stripes is necessary to 
stabilize this charge inhomogeneity. A possible mecha- 
nism, not included explicitly m our model, an in-plane 
long range Coulomb repulsion,^ could help to reduce the 
requirement of an external agent in order to make stripes 
stable. A long-range potential cannot be implemented in 
the world-line QMC used here (it cannot be implemented 
in DMRG calculations either) but it could be included in 
the SEHS scheme as we noticed in the previous section. 

IV. CHARGE AND SPIN ORDERING. 

There is now mounting experimental evidence that 
spin ordering, signalled by incommensurate (IC) peaks 
in the magnetic structure pfacJxu:, follows charge ordering 
as temperature decreases oEj'Ej That is, charge degrees 
of freedom are the driving force of stripe formation. -Sim- 
ilar behaviocjOccurs in stripe formation in nickelatedj and 
manganites.Q In the model here considered, in spite of a 
sizable charge inhomogeneity at high temperatures, IC 
magnetic structure appears at much lower temperatures 
in agreement with experiment. This is a nontrivial result, 
since some theoretical models of stripe forma±km do not 
properly account for this ordering sequence. EJ'til 
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FIG. 3. Phase diagram in the temperature-on-site poten- 
tial plane obtained for the 8x8 cluster, x — 1/8. Tco and 
Tso correspond to the crossovers in the charge and spin sec- 
tors discussed in the text. The point indicated with a star 
was obtained on the 12 x 12 cluster, x = 0.083). Diamonds 
indicate the temperature at which the inverse kinetic energy 
along to the stripes are minimum. 

As an introduction to this section, we reproduce the 
phase diagram inj-the temperature-eg plane obtained in 
our previous worktj for SC stripes (Fig. ]^ . By measuring 



charge and spin static structure factors (C(k) and S'(k), 
respectively), two crossovers have been determined. At 
high temperature there is a crossover in the charge sector 
determined by a change in the peak of C(k) from (0,0) 
to (2(5, 0) {S — 2xtt). At a much lower temperature, there 
is a crossover in the spin sector signalled by the peak of 
S'(k) splitting from (7r,7r) to (7r±(5, tt), which correspond 
to the peaks seen in neutron scattering experiments. Be- 
tween these two crossovers, typical temperature or energy 
scales related to transport properties are located, as we 
will discuss in Section ^. These temperature scales have 
been adopted Ijo experimentally determine the onset of 
the pseudogap.Q 

The change of the peak in C(k) in our model is just 
a consequence of stripe formation: (2(5, 0) corresponds to 
the main Fourier component of a pure SC stripe (away 
from (0,0)). The change of the peak in S^) corre- 
sponds to the formation of anti-phase domains,cl a feature 
which is reproduced in our model.tj In the following sub- 
sections we will more systematically study the effect of 
the anisotropy in the exchange term of the Hamiltonian 
and the doping dependence of the results for SC stripes. 
In the last subsection, we report similar studies for BC 
stripes. 

A. Effect of the XY term. 

Although an antiferromagnetic (AF) correlation in 
the direction perpendicular to the Cu-0 planes would 
imply as a second order process an effective exchange 
anisotropy in the plane, it is clear that the most interest- 
ing model to study is the one that involves an isotropic 
(or nearly isotropic) Heisenberg term. Unfortunately, 
a more isotropic exchange implies stronger minus sign 
problem and hence, we are limited to consider 7 < 0.5. 
However, as we show below, one can get an insight on how 
and to what extent a more isotropic coupling changes the 
results obtained on the Ising limit. 

Let us consider the hole-hole and spin-spin correlations 
defined by 

C(i,j) = {n,n;) (3) 

and 

^(i,j) = {StSD (4) 

respectively, ni is the hole occupation number at site i; 
i = {x,y) and j = {x' ,y'). In particular, if a SC stripe 
along y is located at a; = a:;^ , we are interested in the cor- 
relations along the stripe, where i = {xs,Q), j — {xs,r), 
and in correlations across the stripe, corresponding to 
i — {xs — a,y), j = (xs + a,y) (the inset of Fig. |^ shows 
the case of a = 1). 

In Fig. |, the spin-spin correlation across a vertical 
stripe (indicated with filled circles in the inset) is shown 
for the 8x8 cluster, x = 1/8, J — 0.7, Cs = 1.5, and 
7 = 0, 0.25, and 0.5. These correlations show the change 
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FIG. 4. Spin-spin correlation across the stripe (shown in 
the inset) on the 8x8 cluster, x = 1/8, J = 0.7, Bs = 1.5, 
and several values of the XY term 7. 

of sign corresponding of the spin domains going from in- 
phase to anti-phase at the magnetic crossover tempera- 
ture. Their amplitude decrease as expected with increas- 
ing 7 but the crossover temperature remains roughly un- 
affected. In this, and in in similar figures, typical error 
bars are about the size of the largest symbols used. 
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FIG. 5. Hole-hole correlations along the stripe on the 8x8 
cluster, X = 1/8, J — 0.7, es = 1.5, and several values of the 
XY term 7. Results obtained by exact diagonalization on the 
quarter filled 8-site ring at zero temperature are also shown 
for comparison. 

Another important result of Ref. |l^, viz. the metallic 
behavior of the hole-hole correlations along a stripe is also 
quite robust against the increase of the XY term of the 
exchange interaction. In Fig. ^ we show these correla- 



tions for the 8x8 cluster, x = 1/8, J = 0.7, = 1.5, and 
7 = 0, 0.25, and 0.5 extrapolated to zero temperature. It 
can be seen that their behavior remains unchanged with 
7 except for a uniform shift, within error bars (of the or- 
der of the symbol size) . This shift corresponds to a larger 
filling of the stripe as 7 increases. 



B. Doping dependence. 

The study of the dependence of the results with doping 
is essential to determine the extent to which the present 
mode is physically realistic. For this study, we performed 
QMC simulations on the 12 x 12 cluster, with 12 holes 
{x = 0.083) and 24 holes {x = 0.167). In the first case, 
two SC stripes six lattice spacings apart were stabilized 
by an on-site potential, and in the second case four SC 
stripes three lattice spacings apart. We have in addi- 
tion results for the 8x8 cluster with eight holes and for 
the 12 X 12 cluster with 18 holes both corresponding to 
X — 0.125. In the latter cluster, care should be taken be- 
cause the three SC stripes present imply that the three 
intervening spin regions cannot be in anti-phase. 




FIG. 6. Spin-spin correlations on the 12 x 12 cluster, 
X = 1/8, J = 0.35, Cs = 2.0, 7 = 0, at several dopings 
indicated on the text. The inset shows the spin-spin correla- 
tion along the stripe vs. distance (open circles) at T = 0.09 
compared with exact results on a quarter filled 12 sites chain 
(filled circles) at the same temperature. 

The spin-spin correlations across the stripe are shown 
in Fig. ^. The amplitude of the magnetic order both 
above and below the crossover temperature decreases 
with doping and so it does the crossover temperature 
itself, which is in agreement with experiments. At x — 
0.125, for the 8x8 cluster, the crossover temperature is 
larger than that for the 12 x 12 one (due to the frustration 
mentioned above), but smaller than the corresponding at 
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X = 0.083 (see Fig. 3). The inset shows S{r) along the 
stripe for the same parameters, at the fowest reachable 
temperature (T = 0.09), and the corresponding results 
for a chain obtained by exact diagonalization at the same 
temperature, showing a remarkable agreement between 
them. 

The hole-hole correlations along a stripe, extrapolated 
to zero temperature, are shown in Fig. ^ for the same pa- 
rameters as in Fig. |[ These correlations also present the 
Akp oscillation of a quarter-filled chain seen in Fig. ^ with 
no signs of hole attraction. This behavior, together with 
the above mentioned similarity of spin-spin correlations 
indicate that SC stripes behave at low temperatures like 
an isolated chain. 
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FIG. 7. Hole-hole correlations along a stripe on the 12 x 12 
cluster, X — 1/8, J = 0.35, = 2.0, 7 = 0. Results ob- 
tained by exact diagonalization on the 12 site chain at zero 
temperature are also shown for comparison. 



C. Bond-centered stripes. 

Bond centered stripes are stabilized by applying Csi = 
— Cs < on equally spaced stripes consisting of two near- 
est neighbor columns and es,; = es > otherwise. In 
this case, the sign problem, which is much worst than for 
SC stripes, imposes us further limitations on the accessi- 
ble parameters and temperatures. Hence, this study was 
done for the 8x8 cluster at x — 1/8 filling and on the 
Ising limit of the Hamiltonian. 

If the stripe legs along the y direction are located at 
X — Xg and Xs + 1, then the correlations across the stripe 
involve the sites i = {xs — 1, ?/), j = {xs + 2,y). The rung 
correlations relate the sites i = {xs,y), j = {xs + 
The correlations along the stripes may correspond to the 
same leg (i — {xs,0), j = {xs,r)) or to opposite legs 
{i={xs,0),i^{xs + l,r)). 

Our first important result is related to the appearance 
of anti-phase spin ordering of the intervening domains 



at a temperature at which the charges are already es- 
sentially moving on the stripes, in complete analogy to 
what happens for SC stripes and in contrast with some 
speculations based on an spin-only modelEj This result 
can be inferred from the spin-spin correlations shown in 
Fig. ^, obtained for J — 0.7, — 1.5, 7 = 0. For 
J = 0.35, Cs < 2, the sign problem inhibits us to reach 
low enough temperatures. The variation of stripe fill- 
ing with temperature, also included in this Figure, shows 
that its zero temperature value is mostly achieved al- 
ready at T ~ i. The spin-spin correlation across the 
stripe (Si) experiences a sign change, in this case from 
AF to ferromagnetic, at T « 0.15i, after it reaches its 
maximum AF value which occurs in turn once spin order 
in the intervening spin ladders (as shown along the legs 
at the maximum distance, 5*2, and on the rungs, S3), is 
well established. In fact, the inter-domain correlation Si 
starts to vary with temperature at a much lower temper- 
ature than the intra-domain ones, S2 and S3. The spin 
correlation on stripe rungs {S4) essentially vanishes be- 
low the crossover temperature. These combined features 
suggest that, at least for the parameters examined, the 
formation of anti-phase domain is a result of a collective 
interplay between spin domains and not a result of local 
correlations on the stripes. 
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FIG. 8. Spin-spin correlations on the 8x8 cluster, a; = 1/8, 
J = 0.7, Es = 1.5, 7 = 0, with bond centered stripes. 
Si,i = 1,... ,4 are defined in the text. The hole density 
on the stripe, < Ustr > is shown with crosses. 

On the other hand, the peak in the charge structure 
factor at (7r/4,0) remains unchanged down from T « 5 
for the same parameters of Fig. ^. At larger temperatures 
C(k) is virtually constant at nonzero momentum. 

Let us examine next hole-hole correlations. In Fig. ^ 
we show these correlations, for the same parameters as 
before, as a function of distance at the lowest temper- 
ature available. We have also included for comparison 
the exact correlations obtained for a 2 x 8 quarter-filled 
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ladder at T = 0, with the same parameters and various 
values of 7. Two features are apparent. First, for 7 = 0, 
correlations on the stripe are quite different from their 
zero temperature, isolated ladder, counterparts. In the 
inset of this Figure, it is shown that the temperature evo- 
lution of the stripe correlations is quite smooth down to 
the lowest temperature reached. This behavior suggests 
that BC stripes do not behave like isolated ladders unless 
there is some other crossover at lower temperatures. Sec- 
ond, notice that in isolated ladders pairing occurs near 
the isotropic limit of the exchange interaction (7 > 0.5) 
of the t — J model. Unfortunately, it is not possible to 
reach low enough temperatures and large values of 7 with 
our QMC algorithm to further investigate these issues. 



0.3 




FIG. 9. Hole-hole correlations on the 8x8 cluster, x = 1/8, 
J = 0.7, 63 = 1.5, 7 = 0. Results obtained by exact diago- 
nalization on the 2x8 cluster at zero temperature are also 
shown for comparison. 

V. CHARGE TRANSPORT. 

Another important issue a theory of stripes should 
address is the behavior of resistivity in the under- 
doped region. Resistivity measurements on high- 
quality ia2_a;Sra;Cu04 (LSCO) and YBa2Cu30y single 
crystals^ show a metallic behavior at moderate temper- 
atures in a wide range of doping. These experiments also 
show a divergent resistivity as T — > which was earlier 
detected in LSCO (Ref. ^2|) and in Bi2Sr2_2;La2;Cu06+5 
(Rcf. |3^) even after suppressing superconductivity with 
magnetic fields. The inverse mobility of carriers, also 
shown in Ref. |3^, has a similar behavior with tempera- 
ture as the resistivity. 

A measure of charge mobility is the kinetic energy, 
which can be easily computed in QMC and it is affected 
by relatively small errors. In addition, the kinetic en- 
ergy, due to the use of PBC, is the most important con- 
tribution to the Drude weight (defined below). Although 
a straightforward comparison with the above mentioned 



experimental results is not possible, one can get some 
interesting qualitative insights in the problem. 

Further information on charge transport properties can 
be gained by looking aUthe "paramagnetic" contribution 
to the Drude weight.o To this end, we computed the 
current-current correlations in imaginary time defined as: 

Caa(T) = Oa(r)ja) = ^{jaAT)ja.,m) (5) 

where the paramagnetic q = (0, 0) current operator along 
direction a (a = x^y) is: 

ja = XI " «^ Xl(^i+a,<TCi,^ - C^o-CZ+Q.^t) (6) 

and j„,j(t) = e''"ja,ie-^", t = lA, < I < M. Z = 
Tr{e~^^} is the partition function. Since conserves 
particle number and spin the calculation of {ja.i{T)ja,m) 
is relatively simple. Since, 

Oa(r)j„) = - Y^T^" 

where o-'aregi'^) ^^"^ regular part of the real part of the 
conductivity along the a-direction, it can be shown that: 

Xa= Oa(T)j„)dT = - / a' {uj)duj (8) 

Jo Jo " 

Then, from the standard f-sum rule one obtains the 
expression for the Drude weig 

— = -X. + Uk^) (9) 

TT N 

More elaborated analysis, such as the study of the de- 
pendence of the conductivity with frequency, imply the 
analytical continuation of imaginary time current-current 
correlations to real frequencies by, for instance, the max- 
imum entropy method. This procedure is no longer in 
principle exact and results would be affected with large 
error bars. 

Let us first examine the inverse kinetic energy per site 
as a function of temperature and for various doping frac- 
tions in the case of site centered stripes. In Fig. |l^ 
we show the results of QMC simulations for the inverse 
kinetic energy per site along and perpendicular to the 
stripes, Hy = {Ky/N)~^ and = [Kx/N)~^ resp., on 
the 12 X 12 cluster, J = 0.35, = 2.0, for the t - 
model, and fiUings x = 0.083, 0.125, and 0.167. Both 
quantities decrease with doping, following the general 
trend as the inverse mobilities in Ref. Most inter- 
estingly, it can be seen that although Ky has a positive 
slope in the whole temperature range, consistent with a 
metallic behavior of the stripes and in agreement with 
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the results obtained in the previous section, Kx changes 
it sign from positive slope at high temperature to nega- 
tive at low temperature indicating charge localization in 
the direction perpendicular to the stripes. 
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FIG. 11. Inverse kinetic energy per site along (filled sym- 
bols) and perpendicular (open symbols) to the stripes on the 
8x8 cluster, J = 0.7, = 2.0, 7 = for site centered and 
bond centered stripes. 



FIG. 10. Inverse kinetic energy per site along (filled sym- 
bols) and perpendicular (open symbols) to the stripes on the 
12 X 12 cluster, J = 0.35, = 2.0, and 7 = 0.0. Different 
dopings are indicated on the plot. Dashed line corresponds 
to results obtained on the 8x8 cluster. 

One can determine a temperature scale where and 
Ky start to separate (within error bars) from each other 
as T is lowered. This temperature scale coincides roughly 
with Tco in Fig- Another scale, at a lower tempera- 
ture, can be determined at the minimum of 

i^x^ i.e. at 

the point where charges start to be localized in the direc- 
tion transversal to the stripes. This scale of temperature, 
determined on the 8x8 cluster as a function of Cs has 
been plotted in Fig. |[ From Fig. |o[ this temperature 
scale decreases with doping, following the general behav- 
ior of the pseudogap. This temperature scale is located 
between Tco and Tso in Fig. ^. 

As expected, the effect of the anisotropy in the ex- 
change term of the Hamiltonian is very small for Kx and 
Ky. For 7 = and 0.5, results practically coincide within 
error bars. On the other hand, it is instructive to com- 
pare the behavior for SC and BC stripes. In Fig. ^ we 
show the inverse kinetic energies obtained for the 8x8 
cluster, J = 0.7, Cg = 2.0, 7 = 6, for both types of 
stripes. These results show a larger mobility both along 
and transversal to the stripes and an absence of localiza- 
tion in the case of BC stripes down to the lowest tem- 
perature reached (T = 0.2t). 

The behavior of the current susceptibility Xa is very 
similar to the kinetic energy per site along the same di- 
rection. Typical results are shown in Fig. |lj obtained on 
the 8 X 8 cluster for x = 0.125, J = 0.35, = 2.0, 7 = 0. 



Two differences between those two quantities are appar- 
ent. First, the separation between the kinetic energies 
parallel and perpendicular to the stripes are larger than 
for the corresponding Xy and Xx- Second, both tem- 
perature scales discussed above are larger for the kinetic 
energy per site than for the current susceptibility. Both 
features, which we have seen in all cases examined, sug- 
gest that the diamagnetic response is more sensitive to 
charge inhomogeneity than the paramagnetic one. Notice 
also that finite size effects for these current susceptibili- 
ties are very small. 
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FIG. 12. Current susceptibility and kinetic energy per site 
along and perpendicular to the stripes on the 8x8 cluster, 
X = 0.125, J = 0.35, Cs = 2.0, 7 = for site centered stripes. 
Xx and Xy on the 12 x 12 cluster and for the same doping are 
indicated with and "x" respectively. Xx ^^nd Xv on the 
12 X 12 cluster for x — 0.083 are also shown. 
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VI. CONCLUSIONS 

We have examined various properties of a two- 
dimensional t — J model in which charge inhomogeneity 
is stabilized by an on-site potential by using two differ- 
ent numerical techniques. Diagonalization in a restricted 
Hilbert space and finite temperature Quantum Monte 
Carlo techniques allow us to work with relatively large 
clusters with fully periodic boundary conditions. Extrap- 
olation to the full Hilbert space and extrapolation to zero 
temperature give quite consistent results. 

In the first place, we obtained that the linear filling 
in BC stripes is larger than in SC stripes. In this latter 
case, definitely, a stabilizing mechanism is needed if this 
model has to reproduce experimental results. 

Then, we showed that the effect of a moderate XY 
term in the Heisenberg interaction does not affect the 
main conclusions obtained in the Ising limit in Rcf. 
for SC stripes. First, the anti-phase domain ordering, 
which generates the IC peaks measured in neutron scat- 
tering experiments occurs at a much lower temperature 
than the formation of charge inhomogeneities and charge 
localization. This magnetic ordering is a consequence of a 
collective interplay between the spin domains and are not 
driven by local correlations located on the stripes. Sec- 
ond, hole-hole correlations indicate a metallic behavior 
of the stripes with no signs of hole attraction. The study 
of the doping dependence in the range 0.083 < x < 0.167 
suggest that these features are characteristic of the whole 
underdoped region. 

In addition, we showed that some of the main conclu- 
sions for SC stripes apply to BC stripes as well, although 
in this case our study was mostly limited to the Ising limit 
of the Hamiltonian. Again, as the temperature decreases 
the sequence of events are: first, the formation of charge 
inhomogeneities, then the establishment of AF order in 
the intervening spin domains and finally the 7r-shifted 
magnetic ordering of the whole system. The hole-hole 
correlations on the stripe show again a metallic behav- 
ior and, as a difference with SC stripes, which are very 
similar to isolated chains, they indicate that, at finite 
temperatures, stripes do not behave quite like isolated 
ladders, unless another crossover occurs at lower temper- 
atures than the ones we can reach. 

Finally, the study of charge transport properties pro- 
vides further support to the physical validity of the model 
defined by Eqs. In particular, in the case of site 

centered stripes, we have shown distinctive behavior of 
the inverse kinetic energy, as a measure of charge mo- 
bility, in the directions parallel and perpendicular to the 
stripes, in this latter case showing a characteristic fea- 
ture of localization. A similar behavior is observed in 
the paramagnetic part of the Drude weight computed 
through the current-current correlations. In the case of 
bond centered stripes, on the other hand, the increase 
of the inverse kinetic energy as T — > has not been 
observed down to the lowest attainable temperature. Al- 



though these results were obtained near the Ising limit 
and for well-defined charge inhomogeneities, we believe 
they could help to experimentally decide between SC and 
BC stripes in underdoped cuprates. 
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